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Abstract. We analyze the star formation properties of 
a sample of 21 shell galaxies and 30 early- type galaxies 
members of interacting pairs, located in low density envi- 
ronments (Longhetti et al. 1998a, 1998b). 

The study is based on new models developed to 
interpret the information coming from 'blue' HS/Fel, 
H+K(CaII) and A4000 line-strength indices proposed by 
Rose (1984; 1985) and Hamilton (1985). 

We find that the last star forming event that occurred 
in the nuclear region of shell galaxies is statistically old 
(from 0.1 up to several Gyr) with respect to the corre- 
sponding one in the sub-sample of pair galaxies (< 0.1 
Gyr or even ongoing star formation). 

If the stellar activity is somehow related to the forma- 
tion of shells, as predicted by several dynamical models 
of galaxy interaction, shells have to be considered long 
lasting structures. 

Since pair members show evidence of very recent star 
formation, we suggest that either large reservoirs of gas 
have to be present to maintain active star formation, if 
these galaxies are on periodic orbits, or most of the pair 
members in the present sample are experiencing unbound 
encounters. 
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Galaxies: evolution; Galaxies: formation; Galaxies: funda- 
mental parameters; Galaxies: interactions; Galaxies: star- 
burst 
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1. Introduction 

This paper is a sequel of a series dedicated to studying 
the star formation history of early type galaxies show- 
ing shells, ripples, fine structures etc. which are taken as 
signatures of past interactions (see Malin & Carter 1983; 
Schweizer 1992; Reduzzi et al. 1996) and complex star for- 
mation histories. 

It is commonly accepted that dynamical interactions 
between galaxies leading to the formation of shell struc- 
tures are in general accompanied by star formation if 
enough gas is at disposal to the interacting system. How- 
ever, there are different views of this subject that emerge 
from dynamical models. In fact, numerical simulations of 
merging/ accretion events with SPH codes, in which shell 
structures can develop, yield contrasting predictions as 
far as the accompanying stellar activity is concerned. The 
models by Kojima & Noguchi (1997) suggest that star for- 
mation in the satellite is turned off well before the shells 
develop. This explains why we see in shell galaxies signa- 
tures of post-star-burst. In contrast, the models by Weil 
& Hernquist (1993) show that star formation occurs (in 
the galaxy center) while the shell structure is forming. In 
both cases the shell structure lasts for time scales shorter 
than as 1 Gyr. 

In alternative to the merging/ accretion hypothesis (cf. 
Barnes 1996 and references therein), dynamical simula- 
tions show that shell structures could also originate from 
"weak interaction" between galaxies (Thomson 1991 and 
reference therein). In this case, shells are expected to last 
much longer than with the previous alternative. 

The first paper of the series (Longhetti et al. 1998a) 
presented the sample of galaxies to be examined, i.e. 21 
shell galaxies and 30 members of isolated interacting pairs 
located in low density environments and derived for their 
nuclear region several line strength indices defined in the 
spectral range 3700 A < A <5700 A. The indices were 
further splitted in two groups: red indices (A > 4200 A) 
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and blue indices (A < 4200 A). The red indices (H/3, Mgl, 
Mg2, etc., 16 indices in total) are according to the defi- 
nition by Worthey (1992) and Worthcy et al. (1994) and 
have been transformed into the Lick- IDS standard system, 
whereas the blue indices are A4000 A, H+K(CaII) and 
tW/Fel according to the definition by Rose (1984, 1985), 
Leonardi & Rose (1996) and Hamilton (1985). Finally, the 
spectra used to calculate the indices had a 2.1 A FWHM 
resolution. 

In the second paper (Longhetti ct al. 1998b) the kine- 
matics of the gaseous and stellar components of the galax- 
ies under investigation were properly studied to correct the 
line-strength indices for velocity dispersion. 

With the aid of this material, a thorough comparison 
of the above indices with those predicted by the fitting 
functions of several authors (Worthey 1992; Gorgas et al. 
1993; Worthey et al. 1994; Idiart & de Freitas Pacheco al. 
1995; Buzzoni et al. 1992, 1994) was made by Longhetti et 
al. (1998c) in a subsequent study. Furthermore, the prop- 
erties of the galaxies in the sample were compared with 
those obtained by Gonzalez (1993), with particular atten- 
tion to the distribution on the H/3 vs. [MgFe] plane. The 
main result was that no difference could be detected in 
this plane between interacting and post-interacting galax- 
ies. In particular, Longhetti et al. (1998c) discussed the 
influence of a secondary burst of star formation induced 
by interaction/accretion events. 

In the present paper we will investigate the capabil- 
ity of the blue indices in unraveling the strength and age 
of the last star forming episode. Indeed determining the 
strength and age of this episode would bear very much on 
our understanding of the formation of shell galaxies and 
the relationship between shell formation and star forma- 
tion events. 

The paper is organized as follows. Section 2 summa- 
rizes the basic properties of the "blue" indices. Section 3 
presents models of Simple Stellar Populations (SSP) and 
describes how we calculate the line-strength indices by 
amalgamating the medium resolution empirical library of 
stellar spectra by Jacoby et al. (1984, J84) with the library 
of low resolution theoretical spectra by Kurucz (1992). In 
addition to this, section 3 compares the line strength in- 
dices directly calculated from spectra with those obtained 
from fitting functions (see also Longhetti et al. 1998c). 
Section 4 presents the line-strength indices for composite 
stellar populations (CSP), i.e. model of galaxies in which 
a young stellar component is added to an old one in dif- 
ferent proportions. This would mimic a recent burst of 
stellar activity taking place in an old galaxy (Leonardi 
& Rose 1996). The effects of emission lines and metallic- 
ity on the intensity of the blue indices are also discussed 
in some detail. The results of these simulations and how 
they compare with other similar models in literature are 
discussed in section 5. Finally, some implications for dy- 
namical models of shell formation and galaxy evolution 
arc discussed in section 6. 



2. The Line Strength Indices 

2.1. H+K(CaII) 

According to its definition (see Longhetti et al. 1998a), the 
index H+K(CaII) is actually a measure of the Balmer He 
line. It is in fact only the blend of He with the H(CaII) 
line that is able to affect the apparent ratio between H 
and K(CaII) lines. 

Rose (1985) has pointed out that stars of type later 
than F5 show a constant ratio between H and K(CaII) 
lines, while hotter stars, with a deeper Balmer series, are 
characterized by a stronger blend H(CaII)+He relative 
to K(CaII) and consequently by a smaller value of the 
H+K(CaII)index, which reaches the minimum value in A- 
type stars characterized by the maximum Balmer absorp- 
tion. 

2.2. HS/Fel 

This index, introduced by Rose (1985), is a measure of the 
Balmer US line, whose central intensity is related to the 
average intensity of two Fel lines. At increasing surface 
temperature of a star, this index decreases because the 
H<5 and Fel lines increase and decrease, respectively. The 
H<5/Fel index reaches its minimum value in AO- type stars 
(whereby the Balmer absorption lines have their maximum 
intensity). It turns out to be a good age indicator (just 
like all the Balmer lines) even though it is also sensitive to 
metallicity via the dependence on the Fel lines. In addition 
to this, like the other Balmer lines indices, it might be 
contaminated by emission lines. 

2.3. A4000 

The spectral break at 4000A is produced by two factors: 
the concentration of Balmer lines toward the asymptotic 
limit of 3650A and the increase in stellar opacity at short- 
ward of 4000A caused by ionized metals. The A4000 in- 
dex decreases with increasing surface temperature, it de- 
pends upon the stellar gravity only for stars cooler than 
F5 (Hamilton 1985) and is also sensitive to the chemical 
composition (Van den Bergh 1963, van den Bergh & Sack- 
mann 1965, Hamilton 1985). This index provides informa- 
tion about the stellar population at the turn off, blending 
metallicity and age effects (Worthey 1992). Furthermore, 
the A4000 index is useful for studies of distant galaxies, 
because thanks to the large band-passes adopted for its 
definition it can be measured with small errors also on 
spectra at low signal to noise. Finally, just like all other 
"blue" indices, cosmological recession shifts this index to- 
ward redder spectral ranges whereby modern CCDs are 
more sensitive. 
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3. Line-Strength Indices of SSPs 

Most sudies of evolutionary population synthesis (EPS) 
are based on low resolution stellar spectra, see for instance 
Bressan et al. (1994, 1996) who made use of the Kurucz 
(1992) library. This fact prevents the direct measure of 
line-strength indices on the integrated spectral energy dis- 
tribution (ISED) of stellar aggregates of any complexity, 
e.g. SSPs and galaxy models. The problem is solved by 
making use of the so-called fitting functions that express 
the intensity of the indices as a function of three basic stel- 
lar parameters, namely effective temperature, gravity and 
metal content (see Worthey 1992, Worthey et al. 1994). 
Knowing these quantitities for each single star in the stel- 
lar mix, the integrated line strength indices can be easily 
calculated both for SSPs and galaxy models in a fully con- 
sistent manner with their associated ISEDs, see Bressan 
et al. (1996) for all details. It goes without saying that this 
is possible only if the fitting functions are known. Unfor- 
tunately, there are indices, for instance the blue ones, for 
which fitting functions arc not available. Therefore one is 
forced to make use of spectra with the required degree of 
resolution. This is the subject of the coming section. 

3.1. Using empirical Medium-Resolution Spectra 

In general, the spectral energy distribution sp\(t, Z), flux 
(in suitable units) as a function of the wavelength, of a 
SSP of any age t and metallicity Z is defined by 

f M max {t) 

s Px (t,Z)= <f>(M)f x (M,t,Z)dM (1) 

JM min 

where the product between the initial mass function (IMF) 
4>(M) and the monochromatic flux, f\(M, t, Z) of a star of 
initial mass M, metallicity Z and age t, is integrated along 
the isochrone associated with the SSP in question. The 
integration is performed from M m i ni the minimum mass 
of objects that can become stars, up to M max (t), which is 
the maximum mass of the stars in the most evolved phase 
still contributing to the total luminosity (see for details 
Bressan et al. 1994, 1996). 

The function f\(M,t,Z) indicates the adopted spec- 
trum of the constituent stars of mass M, age t and metal- 
licity Z, i.e. the adopted spectral library as function of 
effective temperature, gravity and metal content (a triv- 
ial re-normalization of the spectra is required to scale the 
flux emitted by a star according to its total bolometric 
luminosity) . This library has of course its own resolution 
which determines the final result. 

For the purposes of this study, we have adopted the 
empirical library of stellar spectra by J84, which contains 
fluxes in units of [erg cm - 2 s- 1 A~ 1 } with a mean res- 
olution of 4.5 A (FWHM) for stars of different spectral 
type and luminosity class. The spectra cover the range of 
wavelength 3500A - 7500A and refer to stars with solar 
composition. 



To incorporate the J84 library in our population syn- 
thesis algorithm, we proceed as follows: 

(i) We assign to each star of the J84 catalog the gravity 
and effective temperature according to the relationships 
by Garcia- Vargas (1991). 

(ii) As the J84 spectra cover only a limited spectral 
range, we have extended them both at the short and long 
wavelength side by means of the Kurucz (1992) spectra of 
the same composition, gravity and effective temperature. 
The Kurucz spectra have been cut at 3510A and 7427A 
and after suitable re-scaling, patched with the J84 spectra. 

(iii) The spectral ranges used to calculate the scaling 
factors arc [3510A-3550A] and [7390A-7427.2A] at the 
blue and red side of the spectrum, respectively. From a 
theoretical point of view, the blue and red scaling fac- 
tors should be the same. Actually, we have detected dif- 
ferences between the blue and red scaling factors of about 
26%. Various reasons concur to build up this difference: 
(a) inaccuracy of the synthetic spectra; (b) errors in the 
flux calibration by J84; (c) some mismatch in the corre- 
spondence between Kurucz (1992) and J84 spectra via the 
correspondence in gravity and effective temperature. 

The monocromatic flux of individual stars in eq.(l) 
F\(M,T,Z) now reduces to F\(M,T, Z Q ) because the 
J84 library is only for solar composition. However, the 
dependence on the metallicity (chemical composition) is 
not completely wiped out, because it still remains via the 
metal dependence of the isochrones along which relation 
(1) is integrated. Though reduced, the effect of different 
metallicities (chemical compositions) on the final results 
can be noticed. 

In order to compare indices derived from direct mea- 
surements on the spectra with those obtained from fitting 
functions, we must now degrade the spectra of our SSP 
and/or model galaxies (see below) to the resolution of 8.2 
A (FWHM) typical of the Lick-IDS spectra. To this pur- 
pose we convolve our spectra with a Gaussian distribution 
with cr= 2.92A. 

Using the degraded spectra we calculate the line- 
strength indices by means of one of the relations below 

I a = (l- |g-)<5A (2) 

I m = -2Mog{^-) (3) 

as appropriate to whether we are dealing with atomic or 
molecular indices (see Worthey 1992, Worthey et al. 1994 
and Longhetti et al. 1998a for the meaning of the various 
symbols). 

In several cases, the results we have obtained do not 
perfectly agree with the standard Lick standard system, 
but a small off-set is present. Determining this off-set is 
not possible because there is no sample of stars or spectra 
in common between the Lick and the J84 library. 
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Fig. 1. Indices for individual stars: comparison between 
indices obtained from fitting functions (indicated with ff) 
and those derived from direct measurements on the stellar 
spectra of the Jacobi et al. (1984) library (indicated with 
J84). The comparison with Buzzoni et al. (1992, 1994) 
(open squares) and Idiart & de Freitas Pacheco (1995) 
(open triangles) is also shown 

Finally, we have calculated indices for three large grids 
of SSPs, i.e. at varying metallicity (Z=0.004, Z=0.02, 
Z=0.05), each grid containing 30 values of the age in the 
interval 0.05 to 19 Gyr. 

3.2. Comparing stellar indices 

In this section we check whether stellar indices derived di- 
rectly from the stellar spectra of the J84 library (therein- 
after the Jacobi indices) are fully consistent with the same 
indices but obtained from fitting functions. The compari- 
son is made for all the stars of the J84 catalog. 

In Fig. 1 we display the comparison between the Jacobi 
indices and those obtained from the Lick system fitting 
functions (Worthey 1992, Worthey et al. 1994) as, how- 
ever, modified by Longhetti et al. (1998a, c) in the high 
temperature regime. 

Similar comparison has been made using the fitting 
function of Buzzoni et al. (1992; 1994) for Mg2, Fe5270 
and H/3 (open squares) and those of Idiart & de Freitas 
Pacheco (1995) for Mg2, Mgb and H/3. The results are also 
shown in Fig. 1 , where the open squares are for Buzzoni et 
al. (1992, 1994) and open triangles for Idiart & de Freitas 
Pacheco (1995). Considering all the uncertainties affecting 
the whole procedure, the agreement we get is remarkable. 



Table 1. Comparison of red indices derived from different 
methods. The symbols ff and J84 stand for fitting func- 
tions and Jacobi et al. (1994) library. See the text for more 
details. The SSP indices are given at three ages (in Gyr) 
as indicated. 



M 


ean Value 


of [Index 


;ff) - Index(J84)] 


Index 


Stars 


SSP 


SSP 


SSP 




(J84) 


ft > 0.1) 


(t > 1.0) 


(t > 3.0) 


H/3 


0.64 


0.16 


0.02 


-0.08 


Mgl 


-0.01 


0.02 


0.03 


0.03 


Mg2 


-0.02 


-0.01 


0.04 


0.05 


Mgb 


0.24 


0.26 


0.51 


0.58 


MgFe 


-0.24 


0.14 


0.31 


0.36 


CN 


-0.07 


-0.01 


-0.01 


-0.01 


CN2 


-0.05 


0.01 


0.00 


0.00 


Ca4427 


0.07 


0.10 


0.17 


0.20 


G4300 


-0.31 


-0.24 


-0.07 


0.01 


Fe4383 


-0.28 


0.21 


0.53 


0.63 


Ca4455 


-0.40 


0.69 


0.94 


1.01 


Fe4531 


-0.22 


0.03 


0.16 


0.23 


Fe4668 


1.78 


0.13 


0.21 


0.22 


Fe5015 


0.70 


0.43 


0.57 


0.59 


Fe5270 


-0.05 


-0.03 


0.08 


0.10 


Fe5335 


0.63 


0.22 


0.25 


0.27 


Fe5406 


-0.07 


-0.05 


-0.01 


0.01 



3.3. Comparing SSP indices 

In this section firstly we examine the red indices for SSPs 
and compare them with their counterparts based on fitting 
functions. Secondly we look at the blue indices for SSPs. 
In such a case, no comparison with fitting function analogs 
is obviously possible. 

In Fig. 2 we show the comparison of SSP red in- 
dices with solar metallicity (Z=0.02) obtained from the 
two methods (fitting functions and direct measurement). 
Along each curve the symbols change according to the age 
of the SSP: full squares for 10< t <19 Gyr, open circles 
for 1< t <9 Gyr, full circles for 0.2< t <0.9 Gyr and open 
squares for 0.05< t < 0.1. 

The data displayed in Fig. 2 show that both methods 
yield consistent results. Although almost all red indices 
strongly depend on the metallicity, i.e. their values signif- 
icantly change for small variations in the metal content, 
this effect is not shown here as we are more interested now 
in the mutual consistency of the two methods rather than 
the absolute values of the indices as function of age and 
metallicity. 

Particularly remarkable is the agreement for those in- 
dices that are commonly considered as valuable age indi- 
cators, such as H/3 and G4300. 

Specifically, the values of H/3 derived from the two 
methods at ages older than 1 Gyr show an average differ- 
ence of w 1% with a root mean difference of « 12%. For 
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Fig. 2. Indices for SSPs: comparison between indices in 
the Lick system derived from fitting functions (ff) with 
those directly measured on the spectrum based on the 
Jacobi et al. (1984) library (J84). The ages of SSPs are 
schematically indicated by full squares (10< t <19 Gyr), 
open circles (1< t <9 Gyr), full circles (0.2< t <0.9 Gyr) 
and open squares (0.05< t < 0.1 Gyr). The metallicity of 
the SSPs is solar (Z = 0.02) 



G4300 we find that they differ by « 8% with a root mean 
difference w 16% passing from one method to the other. 
Also in this case, the average difference gets smaller for 
ages older than 1 Gyr (mean difference w 2%, root mean 
difference w 6%). 

A systematic comparison between the Jacobi and Lick 
indices is presented in Table 1 for all indices in common. 

Finally, in the three panels of Fig. 3 we plot the blue 
indices A4000, H+K(CaII) and StyFel as a function of 
the age for SSP with different metal content. The effect of 
this is shown by using different symbols. The same data 
are also listed in Table 2 (available in electronic form only) 
for the sake of general use. The following remarks can be 
made: 

(a) The indices H+K(CaII) and FW/Fel at ages older 
than 2-^3 Gyr tend to flatten out or, at ages older than 
about 10 Gyr, even to reverse the trend. They somehow 
loose ability in deriving the age. The same is true for ages 
younger than 0.3 Gyr. Therefore, there is only a limited 
range in age in which these indices vary monotonically 
with time, i.e. 0.3 to about 3 Gyr. However the effect of 
metallicity is not negliblc bringing an uncertainty on the 
age of about Alogt ~ 0.3 — 0.4 



Fig. 3. "Blue" indices of SSPs as a function of the age 
(t). The points refer to the solar metallicity [Z = 0.02), 
whereas the heavy solid line is for Z = 0.05 and the thin 
solid line is for Z = 0.004. Limited to the case of solar 
composition, the age intervals are put into evidence: full 
squares (10< t <19 Gyr), open circles (1< t <9 Gyr), full 
circles (0.2< t <0.9 Gyr), full squares (0.05< t < 0.1 Gyr) 

(b) The index A4000 works much better than the pre- 
vious ones, because it is monotonic with the age up to ages 
older than 10 Gyr. The uncertainty brought about by the 
metal parameter is however still large (Alogt ~ 0.3 — 0.4). 

4. Modeling Narrow Band Indices of CSPs 

If fine structures and/or signatures of interaction in early 
type galaxies imply additional star forming episodes of 
various intensities and ages, how they would reflect into 
the spectro-photometric properties of these systems? Do 
the line strength indices keep memory of this past activity? 
How can we model this complex dynamical picture from 
the point of view of spectro-photometry? 

A viable approximation to the problem is to consider 
a galaxy of a certain mass, which underwent bulk star for- 
mation in the far past, acquired its own pattern of abun- 
dances and ages for the mix of stars, evolved passively ever 
since and at a certain age hosted an additional episode of 
star formation of suitable duration and intensity in which 
a certain amount of gas was turned into stars. The source 
of gas can be either internal, if any is left over and retained 
after the first initial activity, or acquired from outside dur- 
ing the interaction episode. 

If the initial star forming activity took place on a time 
scale much shorter than the Hubble time, the host galaxy 
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Table 2. (available in electronic form only) The blue indices A4000, H+K(CaII) and H<5/Fel as a function of the age 
for SSPs with different chemical composition. 



can be further approximated to a single entity, in which 
stars span a certain range of age and metallicity. In the 
context of the standard galactic wind scenario (Larson 
1974), the bulk star forming activity is terminated within 
the first Gyr and the mean metallicity of the stars (galaxy) 
is from nearly two times solar to a third of solar as the 
galaxy mass decreases from 1 to 0.01 xl0 12 M Q (see for 
instance the complete galaxy models by Bressan et al. 
1994, 1996; Tantalo et al. 1996, 1998). 

The bursting mode idea has long been around starting 
from the pioneer studies by Huchra (1977) and Larson & 
Tinsley (1978) to explain the scatter in broad band colors 
observed in normal and peculiar galaxies of the same mor- 
phological type, till the more specific attempt by Leonardi 
& Rose (1996) who investigated the effects of a burst of 
star formation on the indices H+K(CaII) and YiS/Fel . 

Leonardi & Rose (1996) start from the spectrum of a 
normal galaxy, obtained as the average spectrum of about 
70 elliptical galaxies observed in the core of the Coma 
cluster and add to this template the synthetic spectra of 
young SSPs with solar metallicity and different ages. 

The age of the SSP stands for different burst ages, 
whereas the burst strength is simulated by changing the 
relative proportions in which old and young spectra concur 
to build the final spectrum of the bursting galaxy. 

The indices H+K(CaII) and H<5/Fel are then measured 
on the composite spectrum of the post-star-burst models 
and are compared with the observational data. The loca- 
tion of a sample of E+A galaxies on the H+K(CaII) vs. 
H(5/Fel diagram provides a first estimate of the ages and 
intensity of the star forming activity characterizing E+A 
galaxies. 

In this work we make use of the same technique. The 
main difference with respect to Leonardi & Rose (1996) 
is that also for the normal galaxy we adopt a theoretical 
spectrum of suitable age and metallicity. In all experi- 
ments we are going to present the age of the quiescent 
galaxy is 15 Gyr and the metallicity is solar. Tests how- 
ever have been made at varying the age and metallicity of 
the host galaxy. Specifically, we have also adopted for the 
age the values of 12 and 18 Gyr and for the metallicity the 
values Z = 0.004 (20% of the solar value) and Z = 0.05 
(2.5 times the solar value). As expected the results are 
little affected by the choice for the age of the host galaxy, 
whereas they are more sensitive to the metallicity. 

To simulate the bursts of star formation we have added 
to the template spectrum of the host galaxy, spectra of 
SSPs of seven different ages, i.e. 10, 6, 4, 2, 0.9, 0.1 and 
0.05 Gyr. 



In each simulation, the spectra of the two components 
(host galaxy and SSP) are weighted by a factor expressing 
the percentage of stellar mass, relative to the total mass, 
created during the corresponding star formation episode. 

Table 3 shows the correspondence between the frac- 
tionary mass in the burst and the contribution to the to- 
tal luminosity by the two component. This correspondence 
has been calculated using the mass to light ratio charac- 
terizing each SSP at the chosen age. 

The mass is given by the total mass in living stars plus 
the mass in White Dwarfs and Neutron stars for which we 
have assumed the typical values of 0.6 M and 1.4 M©, re- 
spectively. The dependence of the White Dwarf and Neu- 
tron Stars mass on the progenitor mass is neglected here 
for the sake of simplicity. 

The luminosity factor is taken from the Padua library 
of SSPs (Bcrtelli et al. 1994, Bressan et al. 1994, Tantalo 
et al. 1996). 

The percentages of mass supposedly stored in the 
young SSPs ranges from (only the host galaxy is con- 
sidered) to 100 (only the young population contributes to 
the spectrum). 

For purposes of illustration, we show in Fig. 4 a few 
cases with different burst age and intensity (percentage of 
mass turned into stars). Each panel contains the spectrum 
of the hosting galaxy, the spectrum of the SSP represent- 
ing the burst and the resulting total spectrum. Worthy 
of note is that the contribution from the young compo- 
nent grows at decreasing age of the burst and that even 
very small percentages of a young component may deeply 
alter the final spectrum. This is particularly evident look- 
ing at the spectral region short-ward of 4000A, where, for 
instance, the contribution by a SSP of 0.9 Gyr involving 
only the 5% of the total galaxy mass parallels that of the 
host galaxy, 15 Gyr old and containing the remaining 95% 
of the mass. Furthermore, a stellar population of 0.1 Gyr 
representing only the 5% of the total mass can produce 
2/3 of the total luminosity emitted at A < 4000A. This 
comparison makes evident the usefulness of the spectral 
region short-ward of 4000A to unravel even small traces 
of recent star forming activity. 

5. The blue indices diagnostic 

In this section we compare the results of our simulations 
with the data for real galaxies making use of two diagnos- 
tic planes, i.e. H+K(CaII) vs. RS/Fel and H+K(CaII) vs. 
A4000, displayed in Figs. 5 and 6, respectively. 

The galaxy data are from the catalog of Longhetti et 
al. (1998a) further selected by choosing only those objects 
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Fig. 4. Integrated spectral energy distribution (bold lines) of a mixed population made of an old (15 Gyr; thin solid 
line) and a young SSP (dotted lines) component born in a recent star forming event. Left panels refer to the case of 
a young component of 2 Gyr. The percentages of the mass involved in the star-burst episode are indicated at the top 
left of each sub-panel. Right panels show the same but for a young component with age of 0.1 Gyr. In all cases the 
metallicity is solar 



with index H+K(CaII) measured with less than 0.25 units 
of error. Due to the low S/N in the blue part of the spec- 
trum of some galaxies, only 18 shell galaxies and 16 pair 
members pass the scrutiny. The basic data for this sub- 
sample of objects is presented in Table 4. 

The models (dashed lines) superposed to data in 
Figs. 5 and 6 refer to the solar metallicity and show the 
expected values for the three indices in post-star-burst 
galaxies. All lines start from a common origin which is 
the locus in these diagrams of the host quiescent galaxy 
(15 Gyr old and solar metallicity). The indices of this 
reference galaxy are H+K(CaII)~1.21, Hd/Fel~1.12 and 
A4000~2.32. Each line refers to a specific age of the burst 
event as annotated along the curves. Moving from the 
common origin, along each line what is varying is the per- 
centage of mass turned into stars during the burst event. 
The percentage increases moving from the origin to the 
blue end of each curve. This in fact represents the case of 
a 100% young galaxy with the age indicated. 

In these diagrams, an old quiescent galaxy suffering 
from a burst of star formation would soon leave the orig- 
inal place and move toward the blue corner up to a max- 
imum location dictated by the percentage of newly born 
stars. As soon as the burst is over and the new stellar 



component start aging, the galaxy goes back to it original 
place. Once about 3 Gyr have elapsed, the current po- 
sition is virtually indistinguishable from the original one 
(see also the run of the three indices as a function of the 
age in Fig. 3). 

The remarkable novelty showing up both in Figs. 5 and 
6 is the different distribution of the shell and pair galaxies. 
In fact, while shell galaxies fall in the region covered by 
the theoretical models, pair galaxies have a much broader 
distribution. 

The large error bars affecting the data (in particular 
the index FW/Fel) do not allow us to derive any specific 
conclusion about the age and intensity of the secondary 
episode of star formation for individual galaxies. The diag- 
nostic diagrams have to be interpreted only in a statistical 
sense. 

Recalling that old galaxies suffering from a secondary 
burst of star formation perform the loop in the diagnostic 
diagrams away from the quiescence position on a rather 
short time scale (which in turn depends on the amount of 
mass engaged in the burst) , the different behavior of shell 
and pair galaxies can be explained if the bursting episode 
in the former is statistically older than in the latter. 
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Table 3. Simulations of bursting galaxies. T, M sta r, ^boi and M/L are the age, mass and bolometric magnitude 
and mass to light ratios of the SSPs representing the burst. M star and Mf, ; refer to an ideal SSP whose initial mass 
function in number is normalized to one over the whole mass range. The second half of the Table shows the ratio of the 
luminosity associated to the young SSP to the total luminosity of the system at varying the age and mass percentage 
of the SSP representing the burst. 
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This is also somehow hinted by the large value of 
H+K(CaII), larger than 1.2, that some pair galaxies reach. 
Such high values of H+K(CaII) are expected only if the 
index is contaminated by the He emission line (see below 
for more details). Therefore, the distribution of pair galax- 
ies in the two diagnostic planes is statistically indicative 
of ongoing star formation. In contrast, the distribution 
of shell galaxies seem to suggest that the secondary star 
forming episode is on the average older than 1 Gyr. 

We have checked that the above conclusions arc not 
affected by the age adopted for the host galaxy at quies- 
cence. Indeed similar results hold assuming for this object 
ages of 12 and 18 Gyr. 

5.1. Effect of Line Emission 

In order to support our preliminary suggestion that values 
of H+K(CaII) higher than 1.2 are indicative of contam- 
ination by the He emission line and hence ongoing star 
formation, we have analyzed the spectroscopic atlas by 
Kcnnicutt (1992). 

It contains 55 local galaxies, well distributed among 
the different morphological types. Table 5 reports the re- 
sults of our measurements of the H+K(CaII) and A4000 
indices for this sample. All peculiar and irregular galaxies 
have been excluded from the analysis. 



The resolution of the Kennicutt (1992) spectra is be- 
tween 5A and 8A FWHM, slightly better than the one in 
usage here. While A4000 is not expected to depend on the 
spectral resolution and/or velocity dispersion, H+K(CaII) 
tends to increase as the spectral lines get wider. How- 
ever, the expected increase in H+K(CaII), passing from 
the spectral resolution of Kennicutt (1992) data to ours, 
is very small and within statistical errors. 

From the entries of Table 5, the mean value of the 
index H+K(CaII) for early-type galaxies (E+S0) is 1.19 
(± 0.16). In galaxies of later types the index gets smaller. 
H+K(CaII) indeed decreases in stellar systems richer in 
AO-type stars, i.e. stellar systems with active star for- 
mation. Interestingly enough, the value of H+K(CaII) in 
post-star-burst models with solar (Figs. 5 and 6) and non 
solar metallicities (not shown here for the sake of brevity) 
agrees with the maximum observational value ~ 1.2. 

In contrast, the bottom panels of Figs. 5 and 6 show 
some pair members (RR210b, RR381a, RR387a, RR387b 
and RR397b) characterized by values of H+K(CaII)> 1.3. 
This means that the strength of the blend of H(CaII) with 
the Balmer line He is lower than expected. 

The only possible explanation of this fact is to con- 
sider the contamination of the He absorption line by the 
corresponding emission line. 
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Table 4. Sample of shell and pair galaxies with best measured H+K(CaII), A4000 and H<5/Fel 
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It is worth noticing that none of the above 5 galaxies 
shows a clear H/3 emission line. Although the He emis- 
sion line is expected to have a lower absolute intensity 
with respect to that of the H/3 line, its influence on the 
H+K(CaII) index could be greater than that of H/3 on the 
corresponding absorption feature. Then a small variation 
in the intensity of the He absorption line caused by the 
corresponding emission component could be more easily 
revealed in the H+K(CaII) index, thus making it a good 
indicator of star formation. 

5. 2. Effects of the Metallicity 

We have also calculated two grids of post-star-burst mod- 
els with metallicities different from the solar value. The 
results are shown in Fig. 7. As explained in § 3, models 
of non solar metallicity represent only a first order ap- 
proximation because of the incomplete dependence of our 
simulations on this parameter (only in the isochroncs and 
not in the library of stellar spectra) . The new models are, 
however, much similar to the case of solar metallicity pre- 
sented above, even though they differ in some details. 

(i) The reference galaxy of 15 Gyr age and the 
loci of the bursting galaxies shift to bluer and red- 
der indices passing from solar to Z = 0.004 and 
Z = 0.05, respectively. In more detail, for Z = 
0.004 the reference galaxy has indices H+K(CaII)~1.17, 
H(5/Fel~0.82 and A4000~2.05, whereas for Z = 0.05 it 
has H+K(CaII)~1.23, H5/Fel~1.18 and A4000-2.55. All 
the paths of the bursting galaxies are accordingly shifted. 



Table 5. Data from the atlas of Kcnnicutt (1992) 



Class 


A4000 


H+K(CaII) 


E+S0 


2.08 ±0.23 


1.19 ± 0.16 


Sa 


1.73 ±0.36 


1.06 ±0.07 


Sb 


1.57 ±0.11 


0.91 ±0.10 


Sc 


1.42 ± 0.01 


0.91 ±0.03 



(ii) In the H±K(CaII) vs. RS/Fel plane, different 
metallicities produce almost the same loci for the post- 
star-burst galaxies. 

(iii) A significantly larger variation is seen in the 
H+K(CaII) vs. A4000 plane. 

(iv) As compared to the observational data, the Z — 
0.004 models are too blue, whereas the Z = 0.05 ones are 
too red but still marginally compatible. 

Given the adopted approximation and the uncertainty 
in the data (see the error-bars) we cannot contrive the 
mean metallicity of galaxies from their location in the di- 
agnostic planes H±K(CaII) vs. H<5/Fel and H±K(CaII) 
vs. A4000. Nevertheless, the results shown in Figs. 5, 6 
and 7 hint that the solar metallicity is perhaps suited to 
represent the mean value in early-type galaxies. Similar 
conclusions have been reached by Bressan et al. (1996), 
Grcggio (1996) and Longhetti et al. (1998c) analyzing the 
position of early-type galaxies in the H/3 vs. [MgFe] dia- 
gram. 
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Fig. 7. Effects of metallicity on the diagnostic diagrams H+K(CaII) vs. H<5/Fel and H+K(CaII) vs. A4000. The dotted 
and solid lines are for Z = 0.004 and Z = 0.05, respectively. 
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Fig. 5. H+K(CaII) vs. H<S/FeI diagram for shell galaxies Fig. 6. Same as Fig.5 but for the H+K(CaII) and A4000 

(top panel) and pair members (bottom panel). The dashed indices 
lines represent models of post-star-burst galaxies with Z = 

0.02. Each line is labeled by the burst age. See the text 6 - Summary and concluding remarks 

for more details. In each panel the average error-bars are Adopting the medium resolution spectral library obtained 

shown by J84 and the SSP models by Bressan et al. (1994, 1996) 
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and Tantalo et al. (1996) we have calculated three narrow 
band indices in the spectral range A < 4200A, namely 
H+K(CaII), H<5/Fel and A4000, as a function of age and 
metallicity. 

With the aid of these results, we simulated post-star- 
burst galaxies by means of the simple recipe: an old host 
galaxy (with age of 15 Gyr and solar composition) to 
which at some arbitrary age a burst of star formation is 
added. The intensity of this is measured by the percentage 
of mass with respect to the total which is turned into stars 
by the burst episode. 

We have shown that the line strength indices 
H+K(CaII), FW/Fel and A4000 respond to even small 
traces of past star formation activity. In particular, we 
have discussed the sensitivity of H+K(CaII) to the ongo- 
ing star formation events due to the contamination of the 
B aimer He emission line. 

From the analysis of the location of shell and pair 
galaxies on the diagnostic planes based on the indices 
H+K(CaII), HS/Fel and A4000 we may advance some 
suggestions, whose validity is, however, of mere statisti- 
cal nature owing to the large uncertainty still affecting 
the data. 

1. Shell galaxies and early- type galaxies members of pairs 
in the diagrams show different distributions in the 
A4000 vs. H+K(CaII) and FM/FcI vs. H+K(CaII) 
planes. 

2. In shell galaxies, the age of the last star forming event 
goes from 0.1 to several Gyr and involves different per- 
centages of mass. If the last burst of stellar activity 
that affects the line strength indices, correlates with 
the dynamical mechanism, i.e. merger or weak inter- 
action, forming the shell features, in such a case these 
latter are long lasting phenomena. 

3. The distribution of early-type galaxies members of 
pairs suggests that the vast majority of them contain 
a very young stellar component (between 0.1 and <1 
Gyr), i.e. they suffer from a very recent burst of star 
formation. 

Let us now discuss the above results in the context of 
dynamical models of galaxy interaction, having in mind 
that the formation of shell structures of long duration is 
a sort of constraint hinted by the observational data. 

Among the N-body and/or SPH dynamical models 
dealing with the formation of shells, merger models are in 
general unable to produce shells of long duration. Dupraz 
& Combes (1986) may produce shells that last longer than 
in other merger models (Quinn 1984; Hcrnquist & Quinn 
1989), depending on the initial conditions. However these 
models fail in predicting the radial distribution of shells in 
the best studied shell galaxies. Only the inclusion of the 
dynamical friction (Dupraz & Combes 1987), that slows 
down the merging process between the two galaxies, may 
produce long lasting shells with the correct distribution. 



Nevertheless, these models predict also the presence of a 
double nucleus in the galaxy hosting the shells. Notewor- 
thy, in the present sample of shell galaxies only ESO 240- 
100, out of twenty objects, is characterized by the presence 
of a double nucleus. 

Furthermore, SPH simulations of mergers by Weil & 
Hcrnquist (1993) show that during the shell structure 
phase, a large fraction of the gas content (if present) falls 
soon onto the nucleus whereby star formation is likely to 
occur. In these models shells are as old as the last cen- 
tral star forming event. Similar predictions arc found in 
the models by Kojima and Noguchi (1997), in which the 
time delay between the end of the star-burst and the shell 
formation is very short (< 1 Gyr). This prediction hardly 
matches the occurrence of shell structures and relatively 
old star-bursts indicated by line strength diagnostic. 

Alternatively, dynamical models producing shells via 
weak interaction event (Thomson & Wright 1990, Thom- 
son 1991) predict much longer lifetimes for the shells up 
to w 10 Gyr. This type of model requires, however, the 
presence of a thick disc structure in the host early-type 
galaxy. Furthermore, shells do not develop if the galaxies 
experience multiple interactions. 

Rampazzo et al. (1998b) analyzed the isophotal struc- 
ture of pair galaxies in our sample and found that a large 
fraction of them show a disc-like structure. Therefore, if 
the interaction generating shells is the weak interaction 
mechanism of Thomson & Wright (1990) and Thomson 
(1991), we would expect to see them associated to pair 
galaxies where the disc-like structure is often present. 
Among the pair galaxies of our sample, only a few possess 
shell structures, namely RR 225a, RR 225b and RR 278a 
to which the system NGC 1549+NGC 1553 (the differ- 
ence among the two systemic velocities is AV 31 km s _1 : 
Longhetti et al. 1998b) in the Dorado group can be added. 

Finally, let us comment on the observational hint that 
despite the large uncertainty affecting the indices we have 
derived, many pair galaxies seem to be characterized by 
very young bursts. 

Among others, RR278a is particularly remarkable as it 
has the H+K(CaII) index as high as 1.95 and most likely 
H/3 lines in emission, suggesting ongoing star formation. 
Another interesting case is NGC 1553 in the center of 
which Trinchieri et al. (1997), using narrow band imaging, 
have found Ha emission. Longhetti et al. (1998a,b) did not 
detect H/3 emission probably because it was too faint or 
it was filling the absorption feature. If so, also the value 
of the H+K(CaII) line may be contaminated by He in 
emission. It follows from all this that burst in NGC 1553 
could be as recent as in pair galaxies with shells. 

What the implications are as far as the star formation 
history in pair galaxies is concerned? 

Prior to any other consideration, one has to clarify the 
question whether the pair galaxies in question are bound 
or unbound systems. The criteria used to select pairs of 
galaxies are designed to isolate objects whose interaction 
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potential energy is by far larger than that with any other 
nearby galaxy. Nevertheless, there remains the controver- 
sial problem of the relative percentage of bound vs. un- 
bound pairs (see e.g. Junqueira & de Freitas Pacheco 1994 
and reference therein). 

In recent years, efforts have been made to disentangle 
the problem of the physical reality of pairs by means of 
techniques different from the classical statistical methods. 
For instance, the X-ray diffuse background has been pro- 
posed as a probe of the gravitational potential well and 
physical reality of a pair/group in turn. 

The NGC 2300 group, dominated by the pair 
NGC 2300/NGC 2276 (Mulchaey et al. 1993, 1996) and 
surroundend by X-ray diffuse background could be taken 
as an example of a bound system. In contrast, the absence 
of an X-ray diffuse background in K 416 pair of galaxies 
(Rampazzo et al. 1998a), the analog of NGC2300 as far 
as morphology is concerned, suggests this system is an 
unbound encounter. 

In this context, the fact that the majority of pair mem- 
bers in our sample underwent a recent star formation event 
(in some of them even going on now) suggests that either 
large reservoirs of gas must be present to maintain active 
star formation (if these galaxies are bound systems on pe- 
riodic orbits) or most of the pair galaxies in the sample 
are experiencing unbound encounters and undergo an oc- 
casional (perhaps unique) burst of stellar activity. 
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